Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed February 1, 2008 



(MN IMfeX style file vl.4) 



An extended multi-zone model for the MCG— 6-30-15 
warm absorber 



R. Morales 1 , A.C. Fabian 1 and C.S. Reynolds 2 

1. Institute of Astronomy, Madingley Road, Cambridge CB3 OH A 

2. JILA, University of Colorado, Campus Box 440, Boulder, CO 80309-0440 USA 



O 
O 
O 

(N 
X> 



> 

in 
o 

(N 
O 

o 
o 

S3 

Of 

6 ■ 



February 1, 2008 



X 



ABSTRACT 

The variable warm absorber seen with ASCA in the X-ray spectrum of MCG— 6-30- 
15 shows complex time behaviour in which the optical depth of OVIII anticorrelates 
with the flux whereas that of OVII is unchanging. The explanation in terms of a two 
zone absorber has since been challenged by BeppoSAX observations. These present 
a more complicated behaviour for the OVIII edge. We demonstrate here that the 
presence of a third, intermediate, zone can explain all the observations. In practice, 
warm absorbers are likely to be extended, multi-zone regions of which only part causes 
directly observable absorption edges at any given time. 
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1 INTRODUCTION 

X-ray absorption by partially ionized, optically thin mate- 
rial, along the line of sight of the central engine, the so called 
warm absorber, is one of the prominent fea tures in the X- 
ray spectrum of many AGN (Reynoldsl997). The presence 
of such gas was first postulated in order to explain the un- 
usual form of the X-ray spectrum of QSO MR 2251-178 
(Halpcrnl984). A direct confirmation of the existence of cir- 



cumnuclear ionized matter came from ASCA, which for the 
first time was able to resolve the OVII and OVIII absorption 
edges (at rest energies of 0.74 and 0.87 keV , respectively) 
i n the X-ray spectr a of the Seyfert 1 galaxy MCG— 6-30-15 
( Fabian et al.1994 ). Systematic studies of warm absorbers 
in Seyfert 1 galaxies with ASCA have shown their u biquity, 
being detected in half of the sources (Reynolds 1997). 



Warm absorbers were usually described by single zone, 
photoionization equilibrium models. Under this assumption, 
simple variability patterns were expected: increasing ion- 
ization of the matter when t he so urce brightens. This was 
not found by Otani et al. (1996) during their long- look 
ASCA observation of the nearby (z=0.008) Seyfert 1 galaxy 
MCG— 6-30-15. The source varied with large amplitude on 
short time scales (10 s or so). The depth of the OVII edge, 
on the contrary, stayed almost constant, while that of the 
OVIII edge was anticorrelated with the flux. To explain the 
behaviour of the OVII and OVIII edges, the authors adopted 
a multizone model in which the OVII and OVIII edges orig- 
inate from spatially distinct regions. OVII ions in the region 
responsible for the OVII edge, the outer absorber, have a 
long recombination timescale (i.e. weeks or more), whereas 



the OVIII edge arises from more highly ionized material, 
the inner absorber, in which most oxygen is fully stripped. 
The recombination timescale for the OIX ions in the inner 
absorber is of the order of 10 4 s or less. A decrease in the 
primary ionizing flux is then accompanied by the recombi- 
nation of OIX to OVIII, giving the observed variation in the 
OVIII edge depth. 



Orr et al. (1997) raised the possibility of a more com- 
plex warm absorber. During their MCG— 6-30-15 BeppoSAX 
observation the depth of the OVII edge was marginally con- 
sistent with being constant, whereas the optical depth for 
OVIII, t(OV 1 1 1), exhibited significant variability. The au- 
thors claimed that its large value during epoch lfj (1.7±0.5, 
lcr uncertainty) was inconsistent with the values at all other 
epochs. They also pointed out that the value of t(OV III) 
during epoch 5 (a low luminosity state) did not show the 
expected anticorrelaction wit h the ionizing flux . 

In both Otani et al. ( |l996| ) and Orr et al.( [l997| ), 
t(OVIII) was plotted versus count rate. In order to com- 
pare ASCA and BeppoSAX observations, these two obser- 
vational results for r (OVII I) have been plotted versus lu- 
minosity]^] in 0.1-10 keVin figure 1 (instead of versus count 
rate). The conversion from count rate to luminosity is dif- 



* The epochs in the BeppoSAX observation are chronologically 
numerated (i.e. number 1 corresponds to the first epoch of the 
observation, etc.). 

t The conversion from count rate to luminosity has been ob- 



tained using Hr. = 50 km s Mdc " 



_0q_ 



and PIMMS 



(httpj / /heasarc. gsfc.nasa.gov/Tools/w3pimms. html ) 
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Luminosity (10 43 ergs" 1 ) in 0.1-10 KeV 

Figure 1. Comparison ASCA (square) - BeppoSAX (circle) ob- 
servational results for the optical depth of OVIII, t(OV 1 1 1), ver- 
sus luminosity in 0.1-10 keV . Note the only discrepant point in 
the BeppoSAX observation labelled as 1 and the very good agree- 
ment for point 5. 



ferent for each apparatus and therefore, the comparison can 
not be made using count rate, but luminosity. The only dis- 
crepancy between both sets of data is the result for epoch 1 
(labelled 1 in figure 1)^} Note also that there is no disagree- 
ment at all with any other point (Orr et al. mentioned epoch 
5 as problematic, but in this plot it appears consistent with 
the rest of the Otani et al. data). 

In this paper we present a simple photoionization model 
that accounts for the experimental results of both ASCA and 
BeppoSAX observations. Section 2 describes the code used 
to model the inner warm absorber. The application to the 
MCG— 6-30-15 warm absorber is presented in Section 3. The 
extension of the multi-zone model for the warm absorber (i. 
e. the third component) is addressed in Section 4. Finally 
our conclusions are discussed in Section 5. 



2 TIME DEPENDENT PHOTOIONIZATION 
CODE 

In order to reproduce both ASCA and BeppoSAX observa- 
tions, the state of the inner absorber has been modelled 
using a time d ependent photoionization code for oxygen 



(R,eynoldsl996). This code treats the material in the inner 
absorber as containing only ions of oxygen in an otherwise 
completely ionized hydrogen plasma. A given oxygen ion is 
assumed to be in one of the 9 states corresponding to its 
ionization level. Different excitation levels of a given ioniza- 
tion state are not trea ted. The total abundance of oxygen 
is fixed at 7.41 x 10 _4 (Grevessel989) relative to hydrogen. 



t This epoch was at the very beginning of the BeppoSAX ob- 
servation. In a private communication, A. Orr notes that this 
high value for the optical depth for OVIII could not be due to 
instrumental effects, since no similar behaviour, neither at the 
beginning nor during the observation has been observed in any 
other BeppoSAX target yet. 



Both ions and electrons are assumed to be in local thermal 
equilibrium (LTE) with a common temperature T and the 
plasma is assumed to be strictly optically-thin at all frequen- 
cies (i.e. radiative transfer processes are not considered). 

A point source of ionizing radiation is placed at a dis- 
tance R from this material with a frequency dependent lumi- 
nosity L v (total luminosity L). In the calculations presented 
here, the ionizing continuum is taken to be a power-law with 
a photon index F = 2 extending from v min to v m ax- i.e. we 
take 



Au 



(1) 



where A = Yn.{y max /v m in)- The upper and lower cutoff 
frequencies are chosen such that hv max = 100 keV , and 
A = 10, corresponding to hv m in ~ 4.5 eV . 

The physical processes included in the code are: pho- 
toionization, 

Auger ionization, collisional ionization, radiative recombina- 
tion, dielectronic recombination, bremsstrahlung, Compton 
scattering and collisionally-excited OA1035 resonance line 
cooling. Given these processes, the ionization and thermal 
structure of the plasma is evolved from an initial state as- 
sumed to have a temperature of T = 10 5 K , and equal 
ionization fractions in each state. 

The ionization structure of the oxygen is governed by 9 
equations of the form: 



drii 
~~dt 



(2) 



prod. 



where m is the number density of state-i, S stands for 
the ionization/recombination rate per unit volume, the first 
summation on the right hand side (RHS) is over all ion- 
ization and recombination mechanisms that produce state-i 
and the second summation is over all such mechanisms that 
destroy state-i. 

The thermal evolution of the plasma is determined by 
the local heating/cooling rate and the macroscopic con- 
straint. The net heating rate per unit volume is given by: 



dQ 
dt 



(3) 



where A stands for the heating/cooling rate per unit vol- 
ume, the first summation on the RHS is over all heating 
processes and the second term is over all cooling processes. 
For a constant density plasma, we have du = dQ where u is 
the internal energy per unit volume and is given by: 



M = -n e fc_gT. 



(4) 



Thus, the rate of change of temperature is related to the 
heating/cooling rate via: 



dT 
dt 



dQ 



3n e fcs dt 



(•») 



The state of the plasma is evolved using the ionization 
equations represented by (2) and the energy equation (3). A 
simple step-by-step integration of these differential equation 
is performed. The time step for the integration process is 
allowed to dynamically change and is set to be 0.1 x t s h 
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Figure 2. Comparison of the model (solid curve) with CLOUDY 
(dashed line) of a pure hydrogen/oxygen plasma. 




f/T (erg cm s K ) 

Figure 3. Comparison of the model (solid curve) with CLOUDY 
(dashed line). 



where t 3 h is the shortest relevant timescalefj. In the case of 
a constant ionizing luminosity, it is found that the system 
always evolves to an equilibrium state, as expected. 

The model has been compare d against the p hotoion- 



ization code CLOUDY ( Ferlandl996) (version 9004, ( Ferland 
et al.l9^8|)) in the case of a plasma in photoionization equi- 
librium. Figures 2 and 3 report this comparison. There is a 
qualitative agreement between the model (solid line) and 
CLOUDY (dashed line) in figure 2 for the case of a pure hy- 
drogen/oxygen plasma. The main goal of this paper is to 
reproduce the time variability of the oxygen edges, and, this 
figure shows how our model mirrors quite well the behav- 
ior of this element. CLOUDY includes many more physical 
processes and realistic elemental abundances and, as seen in 
figure 3, the agreement is not so good when more elements 
are considered^. We have also investigated the reaction of 
a constant density plasma to an inverted top hat function 
for the light-curve (i.e. initially T = 10 5 K , equal ioniza- 
tion fractions in each state, the luminosity is L — 3 x 10 43 
erg s _1 and held constant. The system is allowed to achieve 
an equilibrium at an ionization parameter £ = — 75 
erg cm s _ 1 , where n is the density of the warm plasma 
and R is the distance of the slab of the plasma from the 



3 At each step the ionization, collisional and recombination 
timescales of each oxygen ionization state is calculated. The short- 
est of these timescales is taken to be t s /,. 

^ Realistic elemental abundances can be important both for 
the ionization structure and for the heating/cooling rate of the 
plasma, and since features of the stability curve result from peaks 
in the heating and cooling functions (which are in turn associated 
with particular ionic species), the comparison presented in figure 
2 is expected to present a much better agreement. The considera- 
tion of realistic abundances is beyond the scope of this paper and 
is currently being studied. 



ionizing source of radiation with isotropic ionizing luminos- 
ity L. The luminosity is then halved for a short period, At, 
and then back to its initial value) . During these changes, 
the state of the plasma is recorded as a function of time. 
This choice of £ places the inner warm absorber a region 
where fos oc L _1 (i.e. in this regime fos is expected to be 
inversely proportional to £ or equivalently to L) . The inves- 
tigation has been carried out for different values of n and R 
keeping n > 2 x 10 7 cm " 3 and R < 1. 4 x 10 17 cm (i.e. the 



constraints given by Otani et al.(1996) for the inner warm 



absorber). Fig. 4 shows an example of the results obtained. 
As expected, fos increases when L is halved, showing a more 
noticeable increment when L is halved over a longer period. 



APPLICATION TO THE MCG 
WARM ABSORBER 



6-30-15 



To reproduce the observational data using our model, it is 
required to obtain the value of the optical depth for OVIII 
once the ionization fraction fos has been calculated by the 
code. Considering the cross section for OVIII, as, as con- 
stant along the line-of-sight and using the fact that the mate- 
rial contains only ions of oxygen in an otherwise completely 
ionized plasma, the optical depth for OVIII, t(OVIII), can 
be written as: 

a s fosARn e 



t(OVIII) 



(6) 



J2 foi x (« - 1) 



where tj% = 10" 19 cm 2 ( |Osterbrockl989| ), AR is the line-of- 
sight distance through the ionized plasma, n e is the electron 
den sity, Abun is t he oxygen abundance relative to hydro- 
gen ( 3reyessel989 ) and foi is the ionization fraction for i 
ionization species for oxygen. 



© 0000 RAS, MNRAS 000, 000-000 



4 



Luminosity 



3x10™ ergs 



Ionization fraction 



°10 4 



2x10* 
Log time (seconds) 



5xl(T 



Figure 4. Time behaviour of the ionization fraction for OVIII, 
/081 f° r different light-curves. The plasma has ionization parame- 
ter £ = 75 erg cm s — 1 for luminosity L = 3 X 10 43 erg s — 1 , elec- 
tron density n e = 3 X 10 7 cm ~ 3 and it is situated at 1.15 X 10 17 
cm from the ionizing source. The luminosity L is halved over an 
interval At equal to: 450(dashed line), 1350(solid line), 4500(dot 
line) and 13500(dot-dot-dot-dashed line) seconds. The corre- 
sponding ionization fraction, /osi is plotted below. 



The light-curve used to rep roduc ed the ASCA obser- 
vation is given by Otani et al. (1996). For the BeppoSAX 
observation, the light-curve has been defined as a step func- 
tion with a constant value for the lumin osity L over the 
different time periods given by Orr et al. (1997). The con- 
stant value for L over each period is chosen to be that given 



in Orr et al. (1997, Fig. 3) plus a period of 2 x 10 s previous 
to epoch 1 with L — 1.8 x 10 43 erg s _1 (a standard value 
for the luminosity during the observation) . Different param- 
eters for the warm absorber have be en in vestigated using 
the constraints given by Otani et al. (1996). With the only 
exception of point 1 (see figure 1), a general good agreement 
for all other experimental points is found. 



Following a suggestion by Orr et al. (1997), we have 



modified the light-curve for the BeppoSAX observation in- 
cluding a previous epoch to it for which the luminosity 
is much lower (the range of values used is 10 40 — 10 43 
erg s _1 )P. The reason for using this very low value for the 
luminosity is that in a regime for which fos ex L~ , the 
highest values fos are expected for low values of L. How- 
ever, the suggested explanation does not reproduce the high 
value of t(OVIII) for epoch 1 in the BeppoSAX observa- 
tion. Even for those values of £ that give a maximum for 
fos (i.e. £ ~ 50 erg cm s^ 1 ), the ionization fraction for 
O(VIII) is still too low to account for the high r(OVIII) 
value at epoch 1. 



II An example of a Seyfert 1 galaxy exhibiting an unusual low 
flux state (a decrease of more than one orde r of magnitude in 
luminosity) is presented in Uttley et al. (11999). 



4 A THIRD WARM ABSORBER 
COMPONENT 

The model we propose to explain both ASCA and Bep- 
poSAX observations incorporates a new zone for the warm 
absorber. Let warm absorber 1 = WA1 be the inner warm 
absorber which parameters are: R < 1.4 X 10 17 cm, n > 
2 x 10 7 cm " 3 and AR ~ 10 14 cm . The outer warm ab- 
sorber will be warm absorber 2 = WA2 with parameters 
R > 3 x 10 18 cm , n < 2 x 10 5 cm ~ 3 and AR ~ 10 14 cm . 
In our model a warm absorber 3 = WA3 is situated between 
WA1 and WA2. The WA3 radius and density will have values 
between those of WA1 and WA2. Therefore, while WA1 and 
WA2 respond on timescales of hours and weeks respectively, 
WA3 is expected to respond to variations in the ionizing 
flux on timescales of days. WA3 would have an ionization 
parameter £ of the order of 500 erg cm s _1 for L — 3 x 10 43 
erg s _1 . Hence the ionization fraction for this value of L 
is too low to be detected. Only when the luminosity of the 
source is sufficiently low (i.e. £ ~ 50 erg cm s ~ ), WA3 re- 
veals its presence by contributing to the total optical depth 
for OVIII. The medium between WA1 and WA2 would be 
constituted by a continuum of warm absorbers: clouds with 
different densities situated at different radii. Only some of 
them happen to be at the radius and have densities that 
efficiently absorbs the OVIII edge ene rgy ( i.e. most of them 
are undetectable). In Baldwin et al. (1995) a model for the 
BLR is presented, in which individual BLR clouds can be 
thought of as machines for reprocessing radiation. As long 
as there are enough clouds at the correct radius and with the 
correct gas density to efficiently form a given line, the line 
will be formed with a relative strength which turns out to 
be very similar to the one observed. Similarly in our model, 
only WA1 and WA2 are detectable for the ordinary values 
of the ionizing flux. Only for the case of a state of low lumi- 
nosity, WA3 will be unmasked. Other zones, as yet unseen, 
may be present. 

Assuming then the presence of WA3 and also an epoch 
of low luminosity previous to the BeppoSAX observation, 
the expected WA3 behaviour would be: 

i) when L « 0.4 x 10 43 erg s _1 (previous to epoch 1), 
the ionization parameter £ ~ 50 erg cm s _1 , giving a high 
value for fos- This period lasts for approximately 10 5 s , so 
the plasma has time to recombine. 

ii) when L ~ (1,5) x 10 43 erg s _1 (i.e. during the ob- 
servation), then £ ~ 150,750 erg cm s~ . For these high 
values of f, oxygen is practically fully stripped and there- 
fore there is a very small contribution to the optical depth 
for OVIII. 

The range of parameters investigated for WA3 is R — 
(2, 8) x 10 17 cm, n = 5 x 10 5 , 10 7 cm ~ 3 and AR in the inter- 
val that gives a column density for WA3 approximately equal 
to 3 x 10 22 cm ~ 3 . After taking into account the soft X-ray 
absorption due to WA1, the response of WA3 has been cal- 
culated and an example of the results obtained is presented 
in figures 5 and 6, where the general good agreement is also 
extended to point 1. 

WA3 has also been modelled using CLOUDY and we 
found a drop in the transmitted portion of the incident 
continuum at ~ 7.8 keV , of approximately 1% (i.e. unde- 
tectable for current instruments). The coronal lines strength 
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Figure 5. Comparison ASCA data (square) with the model 
computations (star) for the following warm absorber parameters: 
WA1: distance to the ionizing source, R = 1.0 X 10 17 cm , line-of- 
sight distance through the warm absorber, AR = 3.0 X 10 14 cm , 



and electron density n e = 4.0 X 10 cm 
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Figure 6. Comparison BeppoSAX data (square) with the model 
computations (star) for the following warm absorber parameters: 
WA1: distance to the ionizing source, R, = 1.0 X 10 17 cm , line-of- 
sight distance through the warm absorber, AR = 3.0 X 10 14 cm , 



and electron density n e = 4.0 X 10 7 



WA3: R = 4.0 x 10 17 



cm , AR = 7.0 X 10 16 cm , and n e = 5.0 X 10 5 cm " 3 . 



has also been checked using CLOUDYp*]. The ratio of the 
modelled to observed fluxes for WA3 are all below 0.1 x f c , 
where f c is the covering fraction^]. 

Finally, the contribution to r(OVII) from each compo- 
nent for our model has been calculated and, as expected, we 
found that WA2 is the main responsible for the OVII edge 
(its contribution makes 96% of the total value of r(OVII)). 



5 DISCUSSION 

Warm absorbers have been the subject of extensive studies 
during the last decade. Such regions are not spatially re- 
solved, and all the available information about their geome- 
try is obtained from analysis of the variability of the oxygen 
edges. The explanation we offer for the time variability of 
the MCG— 6-30-15 warm absorber during both ASCA and 
BeppoSAX observations does not invoke complex processes, 
but a very simple photoionization model together with the 
presence of a multi-zone warm absorber. This would be con- 
stituted by a continuum of clouds at different radii and dif- 
ferent densities, such that only some of them contribute to 
the total optical depth for OVIII depending on the value of 
the luminosity. 

As a final remark note how in our model WA3 is much 
more volume filling than WA1 (AR/R > 0.1 for WA3 while 
AR/R < 10" 3 for WA1). This suggests that WA3 could be 
part of the inter-cloud medium of WA1. 



See Sect. 3.3 (Ferguson ct al.1997) for a discussion of the state 
of the coronal lines atomic data used in CLOUDY, 
tt Porquet et al. (1999) give some restrictions on the density of 



the WA in order to avoid producing coronal line equivalent widths 
larger than observed. Although WA3 presents no problems at all 
for any f c , WA2 does, unless a low value of f c is considered. This 
possibility is currently being investigated. 
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